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Abstract 

In this paper we give a general solution to the problem of the damped harmonic 
oscillator under the influence of an arbitrary time-dependent external force. We em- 
ploy simple methods accessible for beginners and useful for undergraduate students 
and professors in an introductory course of mechanics. 

One of the most important, and investigated, problems in physics is the harmonic 
oscillator. This fact is due to the relevance of this problem in the classical as well as 
in the quantum contexts and, also, because it is one of the few problems we know how 
to solve exactly. The solution to this problem is given in many textbooks of classical 
mechanics, both at introductory and advanced levels, and is given with two arbitrary 
constants determined according to the initial conditions. The treatment of the damped 
harmonic oscillator with the presence of external forces can also be found in the literature 
and its formal treatment, for any time dependent external force, requires the formalism 
of Green's functions. Generally, undergraduate students at introductory levels are not 
yet acquainted with Green's functions formalism, so this subject is usually omitted at a 
first course of mechanics where it is considered only the case of an external force with 
sinus-like dependence. In this last case it is interesting to call attention to the treatment 
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given in Ref. [T] where the problem is solved without resorting to the standard theory of 
second order differential equation. Also, a similar treatment is given in Ref. [2] for the 
single harmonic oscillator. On the other hand, there are also elaborated treatments for 
this problem, particularly interesting is the Feynman diagrams technique given in Ref. [3] 
and the Poisson bracket formalism presented in Ref. [1]. 

In this work, which is intended for undergraduate students and professors, we give a 
general solution to the problem of the damped harmonic oscillator under the influence 
of an arbitrary time-dependent external force. The main advantages of the procedure 
proposed in this paper in comparison with the ones encountered in standard textbooks 
are: the constants are fixed from the beginning, according to the initial conditions of the 
problem, the mathematical methods employed are accessible for beginners once they are 
not sophisticated, we can consider an external force with any time dependence and not 
only the sinus-like one and the method resembles the factorization employed for solving 
the quantum harmonic oscillator. 

The differential equation which governs the time evolution of the damped harmonic 
oscillator with external time dependent force is 

where b is the damping coefficient, k is the restoring constant force, m is the mass of the 
oscillator and F{t) is the external force. The dynamics of the oscillator is determined by 
solving equation ([T]) with the initial conditions, that is, the initial position and velocity, 
x(to) = and x(to) = "^O; respectively. 

In order to solve for x{t) we first define the constants 27 = 6/m and Uq = k/m and 
the differential operator D = d/dt. In this way we write Eq. ([T]) in the form 



(d^ + 27D + ul) x{t) = Fit) , (2) 



which is equivalent to 



[D-a){D-l3)x{t) = F{t) , (3) 

where a and /3 are the roots of the quadratic equation y"^ — 2'yy + ujq = 0, that is, 
a = —7 + iu, l3 = —7 — iu, where u = J ujq — 'y'^ . 
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Using the fact that 



{D - a) fit) = e^'Die-'^'fit)) , 



(4) 



which is vahd for any complex constant a and any differentiable function f{t), we can 
rewrite Eq. (|3]) as follows 

(D - a) [e^'D{e-^'x{t))] = F{t) , (5) 

where we have made a = P and f{t) = x{t). 

Now we use property (jl]) again, but with a = a and f{t) = e^*D^e~^*a;(t)^ , so Eq. 
becomes 

e°*D[e-"* e^*D(e-^*x(t))] = ^ [e-("-^)*D(e-^*x(t))] = e-"*F(t). (6) 
Integrating in the time variable and using the initial conditions we have 



to 



z-'^'' F{t')dt' 



(7) 



where we have performed some simple manipulations on the second term on the left hand 
side of the above equation. 
From Eq. ([7j) we have 



-at' 



to 



F{t')dt' 



(8) 



Performing a second integration and arranging terms we obtain, 



xit) = xoe^^'-'^^ ^V^J^fait-to) _ ^Pit-to)\ ^ 

a — fj ^ ' Jto Jto 



t rt 



"'''F{t')dt' e("-'')*"rft" . (9) 



In order to have just one integration involving the external force F{t), let us manipulate 
the last term of Eq. ([9]). By using the standard expression for integration by parts 



with 



'u{t")dV{t") = U{t)V{t) - U{to)V{to) - fv{t")dU{t") 

to Jto 



U{t") = t e-'^'' F{t')dt' ^ dU{t") = e-'''" F{t")df 

Jto 



(10) 
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dV{t") = e^'^-P^'"dt" V{t") 



a — (3 



(11) 



we can write 

ft rt'" 



/ / e-'^' F{t')dt' e^^~^^' dt" = / e-'^' F{t')dt' 

Jtn Jtn a — 3 Jto 



1 A* 



a — (3 Jto 



f e-^*'F{t')dt' . (12) 

Jto 



Substituting ( IT2|) into Eq. (|9]), using the definitions of a and (3 given after Eq. 
and performing some simple manipulations we have finally 

g-7(i-io) 



X{t) 



UJ 



-It 



Xo 



00 cos(^uj{t — to)) + 7 sm(^uj{t — to)) + ^o sin(^ci;(t — to)) 



2iuj 



iujt I e(T-*^)*'F(t')rft' - e"*^* / e(^+^'^)*'F(t')rft' 



to 



to 



(13) 



Above expression can be written in a more familiar way 



x{t) 



g-7(t-io) 
CO 

+ - / e 

U Jto 



Xo 



oj cos{oj{t — to)) + 7 sin(^ci;(t — to)) + Vq sin(^ci;(t — to^ 



-7(i-i') 



sin(w(t -t'))F(t')rft 



(14) 



in complete agreement with the one encountered in standard textbooks, see for example 
Ref. 0. 



References 

[1] R. Weinstock, "An Unusual Method for Solving the Harmonic- Oscillator Equation" , 
Am. J. Phys. 29, 830-831 (1961). 

[2] R.T. Bush, "The simple harmonic oscillator: an alternative solution of the equation 
for damped oscillation". Am. J. Phys. 41, 738-739 (1973). 

[3] A. Thorndike, "Using Feynman diagrams to solve the classical harmonic oscillator" , 
Am. J. Phys. 68, 155159 (2000). 

[4] CP. Farina and M.M. Gandelman, "An algebraic approach for solving mechanical 
problems", Am. J. Phys. 58, 491-495 (1990). 

[5] K.R. Symon, Mechanics, (Addison- Wesley Reading, MA, 1972), 3'''^ ed. 



4 



